Here we model an intervention to improve food environments and address malnutrition among school-aged children in urban Hanoi. The intervention is formulated as a decision model with three options: introduce a garden at the school, introduce the garden and make it a formal part of the education system, do not introduce a garden. Conceptually the food environment is described in the following figure. The circles with the dotted lines are external or underlying factors.

We express the general concept of the food environment is slightly differently when we describe it in the context of an intervention. In this figure we show that both the garden products and the education program of a school garden are expected to influence the food environment and how this is expected to happen. The garden itself may have a direct effect if it produces fruits and vegetables that are consumed by students. The garden may also effect the food environments of children by influencing their knowledge about and desire for healthy foods, including fruits and vegetables. The education that students receive may change their perception, make them willing to pay, and also increase their access to fresh fruits and vegetables. Several other expected effects are also included in the figure.

Conceptual model of how school gardens in urban Hanoi schools could influence the food environments of children.

Garden model

The gardens will be created in existing land that is available on school grounds. In the the STEM (Science Technology Education and Math) alternative, the gardens will be designed for teaching about nutrition and STEM (an approach to learning and development that integrates science, technology, engineering and maths) at primary and secondary schools.

The model was developed across several iterative workshops in July 2023. These included decision definition, conceptual model development of the selected decision and an initial programming session with preliminary model results. The resulting model was further developed through August to November 2023.

Conceptual model of school gardens as an intervention. Should urban Hanoi school boards invest time and money in creating school gardens? Should they invest in formal STEM education as part of these gardens?

No garden

Here we model the baseline of not having a garden. There are some advantages to this including additional play area for kids and the possibilility of other uses for the space such as parking when school is not in session.

School garden

This intervention involves establishing a garden with a dual focus on gardening education. The initiative encompasses the financial aspects of compost, livestock, garden design, equipment, and construction costs. Establishment costs are allocated for compost initiation, worm maintenance, introducing small animals, hiring a planner for garden design, and acquiring necessary equipment, including potential smart systems. The focus extends from enriching soil quality to creating an aesthetically pleasing and functionally equipped garden

School garden designed to contribute to teaching Science Technology Education and Maths (STEM)

In this option we integrate a focus on sciences and STEM education into the garden. Integrating STEM education through the provision of specialized teaching equipment and training for teachers. Establishment costs include science-oriented teaching equipment, with a focus on practical tools like microscopes. Teachers need training for gardening and STEM topics spanning 5 to 7 subjects over 12 days. The intervention aims to not only create a sustainable and interactive garden space but also to enhance educational experiences by incorporating science-oriented learning.

School garden implementation and monitoring

We want to test our uncertainty here and get as much additional information as we can to better understand the complexity of effects of school gardens on the food environments children. A test run of the intervention, in the form of an actual school garden on school property, will be carried out by the Center for the Development of Organic Agriculture (CODAS) under the Association of Organic Agriculture of Vietnam. In the 2nd year the garden is expected to start running well. The 3rd year is when the STEM education plan will be fully running.

Urban Hanoi school garden

Here is our simulation of the school garden intervention options.

# Source our model
source("CODAS_Garden_Model.R")

# Ensure consistent results with the random number generator
# not for each 'run' of the MC simulation but for 
# consistency each time we call on the simulation 
set.seed(1234) 

garden_simulation_results <- mcSimulation(
  estimate = estimate_read_csv("inputs_school_garden.csv"),
  model_function = school_garden_function,
  numberOfModelRuns = 1000, #run 1000 times
  functionSyntax = "plainNames"
)

Here is a plot of the Net Present Value (i.e. current value of the future benefits) of the garden decision options over 5 years of the intervention. For public and private schools the STEM costs are considered to be in the same garden space but with the additional costs and benefits of a full STEM education program. All options are compared to the same years of using the land for something that is not related to the garden, i.e. as a playground or for parking. Here we plot the distribution for the decision and frame the projected NPV.

source("functions/plot_distributions.R")
plot_distributions(mcSimulation_object = garden_simulation_results, 
                                    vars = c("decision_garden","decision_STEM", "decision_garden_public_school", "decision_garden_STEM_public_school"),
                                    method = 'smooth_simple_overlay', 
                                    base_size = 7, 
                                    x_axis_name = "Comparative NPV outcomes")

This may be easier to see as a boxplot:

source("functions/plot_distributions.R")
plot_distributions(mcSimulation_object = garden_simulation_results, 
                                    vars = c("decision_garden","decision_STEM", "decision_garden_public_school", "decision_garden_STEM_public_school"),
                                    method = "boxplot", 
                                    base_size = 7, 
                                    x_axis_name = "Comparative NPV outcomes")

Summary of results for the decision

Summary

Here we provide a summary of the garden intervention options. We do this with a summary table of the simulation results. We show the percentage of missing values as well as the mean, median and standard deviation (SD) for each output of our model simulations. We use the gt_plt_summary() from {gtExtras} and with options from {svglite}. The table shows the variable name, the plot overview as well as the number of missing values, the mean, median and the standard deviation of the distribution.

# Subset the outputs from the mcSimulation function (y) to summarize only on the variables that we want.
# names(garden_simulation_results$x)
mcSimulation_summary <- data.frame(garden_simulation_results$x[2:70],
 # names(garden_simulation_results$y)
                                 garden_simulation_results$y[1:9])

gtExtras::gt_plt_summary(mcSimulation_summary) 
mcSimulation_summary
1000 rows x 78 cols
Column Plot Overview Missing Mean Median SD
discount_rate 010 0.0% 5.5 5.4 1.5
CV_value 0.010.54 0.0% 0.3 0.2 0.1
inflation_rate 2.812.1 0.0% 7.5 7.6 1.5
size_of_garden 0144 0.0% 52.1 50.7 27.2
expensive_garden_size 88.296.6 0.0% 92.4 92.4 1.5
cost_increase_expensive_garden_size 0.831.79 0.0% 1.3 1.3 0.1
if_students_like 0.360.95 0.0% 0.7 0.7 0.1
if_parents_like 0.311.00 0.0% 0.7 0.7 0.1
if_community_likes 0.01.0 0.0% 0.5 0.5 0.2
if_effective_manage 0.400.77 0.0% 0.6 0.6 0.1
if_garden_yield_enough 0.090.99 0.0% 0.6 0.6 0.2
if_garden_healthy 0.321.00 0.0% 0.7 0.7 0.1
if_teachers_like 0.000.99 0.0% 0.6 0.6 0.2
if_effective_teaching 0.011.00 0.0% 0.6 0.6 0.2
if_effective_training 0.021.00 0.0% 0.5 0.5 0.2
if_offer_green_space 0.261.00 0.0% 0.7 0.7 0.1
if_reduce_pollution 0.090.69 0.0% 0.3 0.3 0.1
if_biophysical_good 0.010.99 0.0% 0.5 0.5 0.2
equipment_cost 126 0.0% 13.1 13.0 4.4
construction_cost 0194 0.0% 58.2 54.0 37.4
garden_designing_costs 7.917.1 0.0% 12.6 12.5 1.5
teacher_training_cost 129 0.0% 12.6 12.6 4.5
school_board_planning 414 0.0% 8.9 8.9 1.8
teaching_equipment 2.612.6 0.0% 7.5 7.4 1.5
compost_starting 112 0.0% 7.5 7.5 1.6
worm_starting 014 0.0% 6.5 6.5 2.2
livestock_costs 036 0.0% 14.9 14.6 5.8
if_family_pays_establishment 0.030.68 0.0% 0.3 0.3 0.1
establishment_family_portion_paid 0.020.99 0.0% 0.5 0.5 0.2
maintaining_labor 1850 0.0% 32.7 32.6 4.5
teacher_salary_cost 1534 0.0% 25.1 25.1 3.1
teaching_equipment_annual 084 0.0% 15.2 11.5 13.1
teaching_tools 013 0.0% 6.0 6.1 2.5
seed_costs 025 0.0% 12.7 12.8 4.4
fertilizer 0.72.6 0.0% 1.5 1.5 0.3
plant_protection 0.86.2 0.0% 3.5 3.4 0.9
livestock_maint 015 0.0% 6.0 6.0 2.5
annual_teacher_training 0643 0.0% 89.8 63.1 87.9
if_school_has_canteen 0.020.66 0.0% 0.4 0.4 0.1
canteen_savings 07 0.0% 3.0 3.1 1.2
sale_of_yield 041 0.0% 20.0 20.0 6.0
extra_cirricular_savings 0151 0.0% 59.0 59.1 24.3
formal_edu_savings 0.14.1 0.0% 2.0 2.0 0.6
formal_edu_savings_STEM 2147 0.0% 59.2 59.1 24.1
outside_investment_value 0.06.6 0.0% 3.0 2.9 1.2
outside_investment_value_STEM 012 0.0% 4.2 4.1 2.1
increased_enrollment_value 010 0.0% 1.7 1.2 1.6
increased_enrollment_value_STEM 0140 0.0% 52.6 50.5 27.6
child_veg_health_care_savings 0.08.7 0.0% 1.7 1.2 1.5
child_veg_school_performance_value 0.000.29 0.0% 0.1 0.1 0.1
child_veg_community_engagement_value 0.000.15 0.0% 0.1 0.1 0.0
garden_mental_health_value 0435 0.0% 76.4 55.8 71.0
child_garden_health_care_savings 0999 0.0% 172.5 122.6 166.3
child_garden_school_performance_value 0260 0.0% 101.4 98.8 49.4
child_garden_community_engagement_value 1.58.7 0.0% 5.0 5.0 1.2
child_STEM_health_care_savings 1120 0.0% 48.9 48.9 17.8
child_STEM_school_performance_value 0175 0.0% 34.9 26.9 30.2
child_STEM_community_engagement_value 0429 0.0% 109.5 94.5 76.3
green_space_value 014 0.0% 5.3 5.1 2.7
reduce_pollution_value 0.33.9 0.0% 2.0 2.0 0.6
school_event_value 01K 0.0% 181.5 136.4 167.6
school_event_freq 030 0.0% 11.1 11.0 5.1
value_of_non_garden_land_use 665 0.0% 35.4 35.7 9.4
if_parking 0.011.00 0.0% 0.4 0.4 0.2
parking_value 0.05.3 0.0% 1.2 1.0 0.9
costs_of_non_garden_land_use 0.16.7 0.0% 2.9 2.9 1.2
land_access 0.041.00 0.0% 0.7 0.7 0.2
suitability_of_land_for_garden 0.361.00 0.0% 0.7 0.7 0.1
beurocratic_barriers 0.070.84 0.0% 0.4 0.4 0.1
NPV_garden -37746K 0.0% 4,956.0 3,104.4 5,584.1
NPV_garden_STEM -2K46K 0.0% 4,582.8 2,707.6 5,599.2
NPV_no_garden 49434 0.0% 250.2 249.3 61.1
NPV_garden_public_school -37730K 0.0% 1,253.9 0.0 3,124.5
NPV_garden_STEM_public_school -2K29K 0.0% 1,027.3 0.0 2,926.3
decision_garden -69546K 0.0% 4,705.8 2,842.8 5,585.8
decision_STEM -2K46K 0.0% 4,332.6 2,500.8 5,601.2
decision_garden_public_school -69530K 0.0% 1,003.7 −217.4 3,127.2
decision_garden_STEM_public_school -2K29K 0.0% 777.1 −229.2 2,929.1
# a summary table with missing, mean, median and sd

Summary of the savings for the passive education garden option

summary(garden_simulation_results$y$decision)
## Length  Class   Mode 
##      0   NULL   NULL

Summary of the savings for the formal STEM education garden option

summary(garden_simulation_results$y$decision_STEM)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
## -1721.1   824.5  2500.8  4332.6  5780.7 45628.3

Summary of costs

summary(garden_simulation_results$y$total_costs)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   86.92  196.88  268.46  406.07  555.20 1964.58
summary(garden_simulation_results$y$total_costs_STEM)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   163.9   311.1   434.6   720.4  1010.0  3654.9

First year

summary(garden_simulation_results$y$Cashflow_garden1)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##  -407.8   105.0   428.3   788.5  1115.4  9300.1
summary(garden_simulation_results$y$Cashflow_garden_STEM1)
##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
## -676.73  -67.99  253.34  606.66  920.52 9018.89

Cashflow of the garden option without formal STEM education

This will be the cost for public and private schools over the intervention.

source("functions/plot_cashflow.R")
plot_cashflow(mcSimulation_object = garden_simulation_results, 
              cashflow_var_name = "Cashflow_garden")

Cashflow of the garden option with formal STEM education

source("functions/plot_cashflow.R")
plot_cashflow(mcSimulation_object = garden_simulation_results, 
              cashflow_var_name = "Cashflow_garden_STEM")

Expected Value of Perfect Information (EVPI)

Here we assess value of information with the multi_EVPI function.

# Subset the outputs from the mcSimulation function (y) by selecting the correct variables be sure to run the multi_EVPI only on the variables that we want. Find them with names(garden_simulation_results$y)
mcSimulation_table <- data.frame(garden_simulation_results$x, 
                                 garden_simulation_results$y[1:9])

Value of information for the garden option (no STEM).

source("functions/multi_EVPI.R")
# first_out_var is the first result variable in the table, "NPV_garden" in our case.
# names(garden_simulation_results$y)
evpi <- multi_EVPI(mc = mcSimulation_table, first_out_var = "NPV_garden")
## [1] "Processing 9 output variables. This can take some time."
## [1] "Output variable 1 (NPV_garden) completed."
## [1] "Output variable 2 (NPV_garden_STEM) completed."
## [1] "Output variable 3 (NPV_no_garden) completed."
## [1] "Output variable 4 (NPV_garden_public_school) completed."
## [1] "Output variable 5 (NPV_garden_STEM_public_school) completed."
## [1] "Output variable 6 (decision_garden) completed."
## [1] "Output variable 7 (decision_STEM) completed."
## [1] "Output variable 8 (decision_garden_public_school) completed."
## [1] "Output variable 9 (decision_garden_STEM_public_school) completed."

Value of information the garden intervention decision.

source("functions/plot_evpi.R")
plot_evpi(evpi, decision_vars = "decision_garden")
## Warning: There are no variables with a positive EVPI. You probably do not need
## a plot for that.

In this case there were no high value variables. We can be confident in supporting the decision, even with our imperfect current understanding.

Value of information for the garden option with formal STEM education.

# using the results of the same multi_EVPI
plot_evpi(evpi, decision_vars = "decision_STEM")

Value of information for the public school garden option with no formal STEM education.

# using the results of the same multi_EVPI
plot_evpi(evpi, decision_vars = "decision_garden_public_school")

Value of information for the public school garden option with formal STEM education.

# using the results of the same multi_EVPI
plot_evpi(evpi, decision_vars = "decision_garden_STEM_public_school")

PLS

We use Projection to Latent Structures model to get some sense of the correlation strength and direction for model variables and our outcome variables.

For passive education garden option

source("functions/pls_model.R")
pls_result <- pls_model(object = garden_simulation_results,
                                resultName = names(garden_simulation_results$y)[1], # the "NPV_garden" 
                                ncomp = 1)
# read in the common input table
input_table <- read.csv("inputs_school_garden.csv")
# source the plot function
source("functions/plot_pls.R")
plot_pls(pls_result, input_table = input_table, threshold = 0.9)

For school garden with formal STEM education

pls_result_STEM <- pls_model(object = garden_simulation_results,
                                resultName = names(garden_simulation_results$y)[2], # the "NPV_garden_STEM" 
                                ncomp = 1)

plot_pls(pls_result_STEM, input_table = input_table, threshold = 0.9)

The full repository can be accessed at https://github.com/CWWhitney/nifam_codas_school_garden with the following QR code.